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ABSTRACT 

Starting from the Lagrangian formulation of the Einstein equations for the vac- 
uum static spherically symmetric metric, we develop a canonical formalism in the 
radial variable r that is time-like inside the Schwarzschild horizon. The Schwarz- 
schild mass turns out to be represented by a canonical function that commutes 
with the r-Hamiltonian. We investigate the Wheeler-DeWitt quantization and 
give the general representation for the solution as superposition of eigenfunctions 
of the mass operator. 
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1. Introduction. 

Recently the dynamics of primordial Schwarzschild black holes has been cast in 
canonical formalism and the quantization procedure has been discussed [1]. A 
complete bibliography that covers the history of the subject is also contained 
there. Indeed, the Hamiltonian formalism is a fundamental key to obtain a quan- 
tum description of a gravitational system and a great deal of work has been 
devoted to the construction of a canonical formalism for the classical black hole 
solutions (see also [2,3]). 

In the present paper we derive the canonical formalism for the vacuum static 
spherically symmetric metric in a simple direct way by foliation in the coordinate 
r. Classically the general vacuum spherically symmetric solution of the Einstein 
equations is locally isometric to the Schwarzschild metric. In order to obtain a 
Hamiltonian description of the Schwarzschild metric we start from the general 
static spherically symmetric line element [4] 

ds 2 = -a(r)dt 2 + N(r)dr 2 + 2B(r)dtdr + 6(r) 2 rf0 2 , (1.1) 

where a, B : N and b are real functions of r and dO, 2 = d9 2 + sin Odtp 2 is the 
metric of the two-sphere. 

Usually, redefining the coordinate time and fixing b = r, (1.1) is cast in the 
form 

ds 2 = -A(r)dt 2 + C(r)dr 2 + r 2 dfi 2 , (1.2) 

where r is now the "area coordinate" since the area of the two-sphere of radius r 
is 47rr 2 . One is then left with two functions A(r) and C(r) that can be determined 
by the Einstein equations. The line element (1.2) is the "standard form" of the 
general static isotropic metric (1.1) [4]. 

The line element (1.1) will be our starting point for a canonical treatment, 
formulated in the coordinate r. We are of course aware that the line element 
(1.1) does not cover the complete spacetime since it describes only a half of the 
Kruskal-Szekeres plane and pure r-coordinate transformations do not lead to 
a complete covering of the Kruskal-Szekeres manifold starting from the metric 
(1.1). In spite of this, the analysis of reparametrizations from this point of view 
may lead to interesting consequences. Indeed, later we will consider a formal 
r-quantization scheme and investigate the ensuing Wheeler - DeWitt (WDW) 
equation [5,6]. 

Since the metric tensor in (1.1,2) does not depend on t, no t-differentiation 
appears in the expression of the minisuperspace action derived from (1.1); starting 
from the Lagrangian we may develop a formal Hamiltonian scheme in the variable 
r and obtain the corresponding r-super Hamiltonian H(a,p a ,b,pi) after having 
introduced the r-conjugate momenta p a and pb- 

Note that there is a range where r is a timelike variable. The signs of iV and 
a are the key. In fact, inside the Schwarzschild horizon of the black hole the area 



coordinate r in (1.2) is a time variable while t is spacelike, and our formalism is 
a true canonical motion in time. In the range where r is timelike, H generates 
the dynamics and plays the role of the usual ADM Hamiltonian; in general the 
r-super Hamiltonian is related to the reparametrizations of the variable r. In 
the metric (1.1) a/|-/V| plays essentially the role of the ADM lapse function with 
respect to the r-slicing [7]. Since we must allow for negative values of N(r), we 
need a slight modification of the ADM formalism, similar to what has been done, 
for instance continuing from a Lorentzian to an Euclidean signature (see e.g. [7]). 

A single Lagrange multiplier imposes the constraint of vanishing of the r-su- 
per Hamiltonian 

H(a,p a ,b,p b ) = Q, (1.3) 

so of course the Hamiltonian "r-dynamics" is generated by a constraint that is 
quadratic in the momenta, as predicted by the ADM canonical formalism. It is 
easy to check that this formalism is equivalent to the Einstein equations for the 
static solution. 

The canonical formalism allows for an interesting algebraic structure of con- 
stants of the motion: in particular we will see that the Schwarzschild mass is 
expressed by a constant canonical quantity, of course gauge invariant. 

The constraint equation H = is independent of r, indeed there has been no 
gauge fixing and r is not determined. The identification of r should be obtained 
by connecting it to the canonical coordinates of the problem (gauge fixing) . This 
procedure can be carried on by the method proposed in [8] for quantum cos- 
mological models. We defer to further study the analysis of gauge fixing and 
quantization in the reduced space. 

We will investigate the quantization of the system by the method of enforc- 
ing the condition H = as an operator condition over wave functions (WDW 
equation). We find the form of the general solution of the equation diagonalizing 
the Schwarzschild mass operator and a commuting operator. The solutions have 
an oscillatory behaviour in the classically allowed regions and an exponential 
behaviour in the classically forbidden ones. 

Thus in this approach the mass plays the role of the quantum number de- 
termining the wave function; in this respect our result is in agreement with the 
conclusions obtained in [1]. 

The outline of the paper is as follows. In the next section we discuss the 
classical r-Lagrangian and r-Hamiltonian formalisms for the metric (1.1). In 
section 3 we integrate the infinitesimal gauge transformations and obtain the 
entire group. We identify the gauge invariant quantities and discuss their algebra. 
Section 4 is devoted to the study of the WDW equation. 

2. Lagrangian formulation. 

Our starting point is the line element (1.1) where the Lagrangian coordinates 
a, 6, B, N are functions of r. As mentioned in the introduction, changes of 



sign in the metric coefficients a and N are allowed (note that the signature is 
Minkowskian over the whole manifold: for instance, if B = 0, aN > 0). r can be 
a timelike coordinate and t spacelike over part of the manifold, so it is a matter of 
preference to define a priori t or r as the timelike variable. Hence, we develop a 
formal canonical structure in r in which the r-super Hamiltonian H is a generator 
of gauge canonical transformations that correspond to reparametrizations of the 
r coordinate in the Lagrangian formulation (and thus in the region where r is 
timelike it generates the dynamics). Hence it seems worthwile to study in detail 
this r-canonical structure. 

Let us consider the line element (1.1) that corresponds essentially to use a 
Gaussian normal system of coordinates with respect to the three-surface (£, 9, <f>), 
i.e. to perform the 3+1 slicing with respect to the r coordinate. As remarked 
in the introduction, looking at (1.1) one realizes that the variable y/\N(r)\ plays 
the role of the r-lapse function in our foliation [7] . The Einstein-Hilbert action 

S =~^ I d^x^R--±- f d s xVhK (2.1) 

16nG Jy 4 8nG J 9Va 

can be cast in the form 

S= r dt f 2 rfrL(a,6,A), (2.2) 



/—la'bb' ab' 2 \ 
L = 2 >/A — - + — + 1 . (2.3) 



where 



(primes denote differentiation with respect to r). In (2.3) we have set AG = 1 
and A is given by 

A(r)=aN + B 2 . (2.4) 

Eq. (2.3) requires that A > and from (2.4) the signature of (1.1) is Minkowskian 
for any value of r. From (2.3) the Einstein equations of motion can be recovered 
considering formally a(r), N(r), B(r) and b(r) as Lagrangian coordinates evolving 
in r. Of course, vA acts as a Lagrange multiplier (and we still have the freedom 
of choosing B(r) or N(r)). From the vacuum Einstein equations derived from 
(2.1), or directly from (2.3), one obtains 

A = k 2 b' 2 , (2.5a) 

a = fc 2 fl-^Y (2.5b) 



where k and M are two integration constants. Since the metric is t-independent, 
we can arbitrarily rescale t in (1.1). This corresponds essentially to fix k in (2.5), 
so we can set k = 1; then the metric coincides with the standard Schwarzschild 



form. M is the Schwarzschild mass. Eqs. (2.5) will be useful for comparison 
with the Hamiltonian formalism that will be developed below. Note that the 
Lagrange multiplier vA can be arbitrarily fixed; furthermore, since A is related 
to N and B by eq. (2.4), also N, or B, can be arbitrarly chosen; these two choices 
correspond to the freedom in the definition of t and r in the line element (1.1). 
For instance, the choice A = 1 corresponds to the area gauge since from (2.5) we 
obtain r = b: 

,1/2 



2M 



ds z = - 1 )dr + N(r)dr z ± 2 



i-(i-^W) 



dt dr + r z dQ, z . 



(2-6) 
The line element (2.6) corresponds to the standard form of the Schwarzschild 
solution for N(r) = (1 — 2M/r) _1 , to the Eddington-Finkelstein metric for N = 
1 + 2M/r and to the line element of ref. [2] choosing N — 1. 

Let us now set up the Hamiltonian formalism in r. We introduce the r— 
conjugate momenta as 

2bb' 
Pa = -/=, (2.7a) 



p 6 = -=(a'b + 2ab'), (2.7b) 



and by the usual Legendre transformation we obtain the density of the action 
(with respect to the coordinate t) 

S = drl-(a'p a + b'p b -ap' a -bp' b )-lH\. (2.8) 

H is the Schwarzschild r-super Hamiltonian 

H = p a (bp b -ap a )-4b 2 , (2.9) 

and 



' = w (2 ' 10) 

has been chosen as Lagrange multiplier. Note that the Legendre transformation 
used to write (2.8) is singular for 6 = 0, but not for a = 0. As a consequence of 
(2.8) we have the constraint 

H = 0. (2.11) 

This constraint expresses the invariance under r-reparametrization and inside 
the region where r is timelike it generates the dynamics. 

Eqs. (2.1) - (2.11) can be easily extended to the Reissner - Nordstrom (RN) 
case, i.e. to a static electrically charged black hole. Let us consider a radial 
electric field whose potential 1-form is 

A = A(r)dt (2.12) 



(this Ansatz was used in [9,10] for the discussion of Euclidean electromagnetic 
black holes). Adding to (2.3) the electromagnetic Lagrangian and using (2.12) 
the Hamiltonian becomes 

#rn = Pa{bp b - ap a ) - Ab 2 + P\ = H + P%, (2.13) 

where Pa is the conjugate momentum to A. Since (2.13) is separable we can solve 
the equation of motion for the electromagnetic field and have Pa = Q where Q 
is the charge of the black hole. Eq. (2.11) becomes 

H RN = H + Q 2 = 0. (2.14) 

The RN case is equivalent to the Schwarzschild case with the constraint (2.14) 
in place of (2.11). 

3. Algebra and Gauge Transformations. 

The gauge transformations of the system are generated by H (i = a, b): 

6qi = a(r)— = a(r)[qi,H] p , (3.1a) 

5pi = -a{r)—— =a(r)[pi,H] (3.1b) 

51 = ^. (3.1c) 

dr 



The action (2.8) is invariant under (3.1) apart from a boundary term that does 
not change the classical equations of motion: 



r-2 



S = I dr — 

dr 



n 



, OH OH , 

v s~i/ri . fin . ' 



dpi dqi 

With a — > l(r)dr eqs. (3.1a,b) are the equations of motion. 



(3.2) 



The system described by the action (2.8) has remarkable algebraic properties. 
Consider the following canonical quantities: 

J = 8b-p a p b , (3.3a) 

/ = b/p a . (3.3b) 

J and / are canonically conjugate gauge invariant quantities (and also obviously 
integrals of the motion): 

[J,I] P = 1, [J,H] p = 0, [l,H] p = 0. (3.4) 

5 



It is also interesting to consider the canonical quantity 

N = bp b - 2ap a . (3.5) 

We have 

N = IJ + 2H/p a , (3.6) 

and the relations 

[N, H] p = -2H, [N, I] p = -I, [N, J] p = J. (3.7) 

N is not gauge invariant, however, in the case Q 2 = 0, i.e. for a Schwarzschild 
metric, it is constant on the constraint H = 0. We shall see that N plays an 
interesting role in the frame of the WDW equation. 

The gauge transformations (3.1) from a gauge l\ to li can be integrated 
explicitly. We have 

a = -Ha 2 +IJa + 4I 2 , (3.8a) 

b = , (3.8b) 

a 

1 

Pa = , (3.8c) 

a 



p b = J a + 8/, (3.8d) 

a = I dr(j,2 — h)- (3. 



From eqs. (3.8) the gauge independent relation follows 



I 2 
a= -2 (4b 2 - Jb-H). (3.9) 



On the constraint H — (Schwarzschild metric) 



a=\I 2 [ 1-^-1. (3.10) 



46 

Therefore from (2.5) 

J = 8M, (3.11) 

where M is the Schwarzschild mass. On the constraint H = —Q 2 the two roots 
of a = in (3.9) correspond to the two horizons of the RN metric. 

6 



It follows that in the case of the Schwarzschild metric I is the momentum 
conjugate to the Schwarzschild mass. This suggests to perform a canonical trans- 
formation to new pair of canonical variables, (J,I;q a ,p a ) where q a = —H/p^. 
This motivates our choice of the eigenfunctions in the discussion of the WDW 
equation. 

4. Quantization. 

The quantization of this apparently simple system exhibits ambiguities that are 
characteristic of the canonical quantization of systems described by general rela- 
tivity [11]. 

A main problem in general is that, in order to set up canonical quantization 
rules, we must know a priori the causal structure of the model representing a 
physical system. To be more specific, we must know which coordinate plays 
the role of time and consequently write down equal time canonical commutation 
relations. 

This is usually an ambiguous procedure. In the classical treatment, the 
identification - if any - of the time variable results from the solution of the 
classical equations of motion and it is not determined a priori. Of course, in some 
cases, as for instance the Friedmann - Robertson - Walker (FRW) model, one 
assumes the signature of the metric (see e.g. [12]). This is because the outcome 
of the equations of motion is anticipated, and a limitation in the signature of 
the metric is consequently assumed. However, strictly speaking, these limitations 
are not always known at the start. This becomes evident whenever the classical 
equations allow for a change in the signature of the metric (see e.g. [7]) or when, 
as in the present case, the presence of a horizon induces a double change of 
signature in the metric. 

In our present case we know from classical solutions that the signature of the 
metric (and the gauge fixing of the coordinate) implies for r a timelike range. It 
is then tempting to explore the implications of a canonical quantization of this 
system imposing equal r commutation relations. This will be carried out in the 
present section. 

We shall impose the constraint (2.11) as an operator condition on the wave 
function. This is the WDW equation. It expresses a necessary condition for 
the wave function, although it does not in general contain all the information 
relevant to the quantum form of the theory. Indeed as it is well known the time 
is not identified, the solution contains both positive and negative frequencies, it 
is a hyperbolic differential operator and thus it does not lead to a well defined 
boundary value problem. It is also plagued by ambiguities since the metric in the 
Hilbert space is not defined. 

We believe that the correct procedure [8] requires identification of the pa- 
rameter r (our internal time) through a gauge fixing condition that defines r in 



terms of the canonical variables and leads to a unitary Hamiltonian in the reduced 
canonical space. When this is possible the quantization of the system is non am- 
biguous and the solutions contain also the information from the constraint. The 
problem of the gauge fixing in the present case will be treated elsewhere; here 
we shall limit ourselves to explore the properties of the solutions of the WDW 
equation. 

The fundamental commutation relations are: 

[a, p a ] = i, (4.1a) 

[b,p b ]=i- (4-lb) 

The operators a, p a , 6, pb have the Schrodinger representation, there being the 
usual ambiguities about the measure to be used. 

We introduce also the mass operator J and its conjugate I according to eqs. 
(3.4). We have 

[I, J] = i. (4.2) 

We remark in particular that J commutes with p a . 

The expression of the WDW Hamiltonian operator is (we consider for sim- 
plicity the case of the Schwarzschild metric, Q — 0) 

H W dw = -apl -bJ + Ab 2 + i\p a . (4.3) 

The term A depends on the ordering and on the representation of p a and p . The 
choice of the covariant Laplace - Beltrami operator [13] leads to A = 1 while the 
symmetric ordering of the operators a,p a and 6, pb leads to A = 1/2. In what 
follows we shall keep A undetermined. 

First of all we determine the eigenfunctions of the commuting operators J 
and p a - We choose the simplest representation, 

P, - -i|. (4-4) 

(q = a, b) . Then the eigenvalue equation for J is 

(8b + d a d b ) Vm = 8MVm, (4.5) 

and the eigenfunctions of J and p a are given by 

-— exp i (pa + p~ 1 /3 M ) , (4.6) 

P 2/T v ' 

where p is the eigenvalue of p a and 

p M =4b(b-2M). (4.7) 



The set (4.6) is orthonormal in — oo < a < +00, —00 < b < +00 with unit 
measure. Let us remark that this approach can be easily adapted to a different 
interval in a, b and to different representations for p a , p^\ the wave function will 
change correspondingly, however the important properties to exploit remain the 
role of J and of the commutation relation [J,p a ] = 0. 

The form of 3m is related to the existence of the horizon at b = 2M for 
positive M. Expressing the solution of the WDW equation 

Hwdw^ =0 (4.8) 

as superposition of ip P M, the general representation of the WDW wave function 
for the Schwarzschild black hole is given by 

»M)=/*^-»/'/* B C( m ) *„(,,»). (4.9) 

C(m) is arbitrary. It is interesting to remark that there is a priori no limitation 
on the sign of the mass m. 

Using a well known representation for the solutions of the Bessel equation 
[14], the representation (4.9) can be cast in the form 

r f 8 \ (A " 1)/2 

*(a,6) = dmC(m)(-^j K 1 . x (2 v / -a/3 m ) (4.10) 

and the solution with fixed mass M is 

/ a \(A-l)/2 

*m = C M\- d f) Ki_a(2 v /Z ^7). (4.11) 

It is natural to assume the form (4.11) of the solution in the regions where a/?M < 
0, namely in the classically forbidden regions a < 0, b > 2M and a > 0, b < 2M, 
where (4.11) is damped exponentially for large b. In the two classically allowed 
regions for the black hole, namely b > 2M, a > and b < 2M, a < 0, the 
behaviour is oscillatory and one should write the appropriate oscillating solutions 
with outgoing or incoming asymptotic conditions. Note that for large b in these 
regions the phase approaches the value of the action evaluated on the classical 
solution for the asymptotically flat spacetime. We are not discussing the joining 
of the wave functions between the different regions as this depends on the choice 
of the ordering and also on the representation assumed for the momenta. 

Suitable superpositions of the kind (4.10) may give wave functions that are 
regular also for 6^0 [10,15]. We note also that the general solution for the Kan- 
towski-Sachs Euclidean wormhole found in [10] corresponds to the solutions of 
the present WDW equation obtained by diagonalizing the operator iV (see (3.5)) 

9 



in place of J . Indeed, using the same choice for ordering (A = 1) and measure in 
superspace as in [10], the differential representation of N is 

N = -t(bd b -2ad a ), (4.12) 

and the solutions of the WDW equation that are eigenfunctions of N with eigen- 
value v are: 



*„(a, b) = \ f 2sE^«A (_ a )W2^ (46v^) • (4.13) 



These solutions are real in the region a < and ortho normal in < b < oo, 
— oo < a < with measure b da db: 

/0 /"O0 

da / dbb^t^ v , =b{v-v'). (4.14) 

Again the phase factor coincides asymptotically with the classical phase factor as 
for (4.11). It is interesting to note that when v = the solution (4.13), namely 

V v=0 =^Ko{4bV=a~), (4.15) 



coincides with (4.11) for M = (and A = 1), as expected since J\l/ M=o( a ,b) = 
iV\I/M=o( a 5 b) = on the constraint shell H = 0. This wave function describes a 
vacuum wormhole in the classically forbidden region. This equivalence supports 
the conjecture [16] that the ultimate remnant in the evaporation process of a 
black hole is a vacuum wormhole. 

5. Conclusions. 

The classical Einstein equations for a static spherically symmetric metric can be 
cast in Hamiltonian form. The starting point is the ADM foliation performed 
along the coordinate r. This is of course a constrained canonical formalism, 
the constraint being that the Hamiltonian vanishes. The Hamiltonian gener- 
ates gauge transformations of the canonical variables that correspond to the 
reparametrization of the coordinate r in the customary formalism of General 
Relativity. 

By a suitable, self - suggesting choice of the Lagrangian multiplier (analo- 
gously to what done in [17,8] for the FRW universe) the Hamiltonian assumes a 
beautiful polynomial form. The infinitesimal gauge transformations can be inte- 
grated, thanks essentially to Einstein and Schwarzschild. This is an interesting 
integrable non linear system. Integrability is due to its simple algebraic structure. 
Indeed, one identifies a pair of conjugate gauge invariant quantities: one of them 
is the Schwarzschild mass. 

10 



Then, the temptation to explore the quantization of this system is big and 
we have carried on the investigation of the WDW equation. In doing this, one 
is comforted by the fact that inside the horizon of a black hole r is a timelike 
variable. 

Note that if we do not fix the coordinate gauge by expressing r in terms 
of the canonical coordinates, this statement is vague: for instance the trivially 
different fixings b = r (area gauge) and b = e r lead to obviously different values 
for the horizon in terms of r. However, this does not matter much: there is a 
region where r is timelike. 

Thus we have studied the WDW equation and give the general representation 
of the solution in terms of superpositions of eigenfunctions of the mass operator. 
It is interesting to observe that there is no reason why the sum should be limited 
to positive eigenvalues of the mass only. 

There is nothing in the form of the WDW equation that reminds us of the 
region in a, b where it is valid, as the WDW equation does not contain r. So we 
may determine the solution in the four regions a^O, b^2M (for positive mass). 
We have not discussed the joining of the solutions between these regions, as the 
result may be affected by the ambiguities in the ordering of the operators and in 
the choice of the measure. 

The solution in the classically forbidden regions can also be cast in a form 
identical to the solution representing a Euclidean wormhole in the Kantowski - 
Sachs spacetime [10]. These solutions are eigenfunctions of a different operator 
N that commutes weakly with the Hamiltonian. In particular the state with 
eigenvalue of N is also eigenstate of the mass with eigenvalue 0. This equivalence 
may support the conjecture [16] that the ultimate remnant in the process of 
evaporation of a black hole is a vacuum wormhole. 

The WDW equation is plagued by the so well known problems. A more 
natural way to investigate the quantum properties of the system seems to be the 
introduction of a gauge fixing of the parameter in the canonical treatment [8] 
that connects r to the canonical variables and leads to a unitary Hamiltonian in 
the reduced canonical space. We defer to a next paper the investigation of this 
method as well as of the connection between the WDW equation and the gauge 
fixed quantization for integrable systems. 
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